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We briefly report on zero-temperature photoassociation of a Bose-Einstein condensate, focusing on 
the many-body rate limit for atom-molecule conversion. An upgraded model that explicitly includes 
spontaneous radiative decay leads to an unanticipated shift in the position of the photoassociation 
resonance, which affects whether the rate (constant) maximizes or saturates, as well as the limiting 
value itself. A simple analytical model agrees with numerical experiments, but only for high density. 
Finally, an explicit comparison with the two-body unitary limit, set by the size of the condensate, 
finds that the many-body rate limit is generally more strict. 
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Introduction. -Quantum degenerate molecules are of in- 
terest for investigations fundamental constants [1] and 
physical chemistry [3J, quantum computing and encryp- 
tion applications [3J, as well as condensed-matter [3] 
and quantum-gravity analogues [5 . For alkali-metal 
atoms quantum degeneracy occurs around 100 nK, 
and is achieved by laser and evaporative cooling atoms 
that begin at an oven temperature of about 600 K. 
Unfortunately, laser cooling depends on the existence 
of a closed few-level system, which does not exist for 
molecules, making quantum degeneracy harder to achieve 
with laser cooling [7]. Coherent association-via either 
laser or magnetic fields-can provide an end run around 
the need for molecular laser cooling: Since it is efficient 
at quantum degenerate phase space densities jS], and 
since it is fundamentally coherent, association creates 
molecular condensate from already-condensed atoms [9]- 
Whereas magnetoassociation is currently in vogue |10j . 
photoassociation is crucial for systems that do not pos- 
sess a Feshbach resonance, such as alkali-earth atoms. 

The idea of coherent photoassociation raises a funda- 
mental question of its own, namely the rate limit for con- 
verting atoms into molecules. Simply put, photoassoci- 
ation depletes the pair wavefunction, and the rate limit 
on atom-molecule conversion is set by the refill time. In 
a many-body model |llj . the refill time is determined by 
the interparticle spacing, and the maximum rate per unit 
density is approximately K = Rj p ~ hp^ 1 ^ 3 /m r , where 
p is the condensate density and m r is the reduced mass 
of a dissociated pair of atoms. In a two-body model [12] . 
the de Broglie wavelength sets the refill time, so that 
K = Rj p r~j T%Kr)/m ri which is also known as the uni- 
tary limit and, at zero-temperature, is set by the size of 
the condensate. The two limits should merge at the onset 
of quantum degeneracy, since pA 3 D ~ 1, but the deeply 
degenerate regime has not been fully investigated. Also, 
both the many-body and unitary rates reach a maximum 
for strong atom-molecule coupling. However, the many- 
body model does not include spontaneous radiative de- 
cay, and this qualitative agreement should be confirmed. 



Experiments above |13) or at |14j the quantum de- 
generate limit are inherently inconclusive, as mentioned 
above, whereas the first experiments with a condensate 
were thwarted by strong dipole forces at strong laser cou- 
pling |15) . Subsequent condensate experiments combined 
photoassociation and the Feshbach resonance, which en- 
hances photoassociation and thereby allowed for observa- 
tion of the rate limit at manageable intensities |16) . but 
observations were again consistent with both the rogue 
and unitary model [17] . Recent theory has shown that, 
up to an interference factor, the many-body result for 
Feshbach-photoassociation agrees reasonably with the re- 
sult for photoassociation alone, but only agrees with the 
unitary limit in a small magnetic-field window near the 
Feshbach resonance, and is otherwise more strict [18] . 

Here we briefly report on the many-body rate limit 
for photoassociation of a Bose-Einstein condensate at 
zero temperature, based on an upgrade to our previous 
model [11 to include spontaneous radiative decay. Nu- 
merical experiments reveal an unanticipated lightshift in 
the position of laser resonance, and the effect of this light- 
shift on the photoassociation rate is examined for high 
and low condensate densities. We use LiNa as an ex- 
ample, but the results apply generally to other species, 
including homonuclear pairs, as well as the Feshbach res- 
onance. 

Model-We focus on the many-body model for heteronu- 
clear photoassociation [19] of a Bose-Einstein condensate 
of two miscible [^H] species of atoms. Photoassociation 
generally occurs on a timescale much shorter than the 
dipole-dipole and trap interactions, and these are there- 
fore negligible. Consider Ni atoms of mass mi and iVjj 
atoms of mass W2 that have each Bose-condensed into 
the zero-momentum state |ai,ci2) = 1 0i ) | O2) , with the 
total number of atoms N — Ni + N2. The photoassoci- 
ation laser then destroys an atom from each species and 
creates a dipolar molecule of mass = mi + m,2 in the 
electronically-excited state |e), with the laser detuning 
6 > indicating an open dissociation channel, as per 
Fig. flTa). The condensates are assumedly miscible, so 
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that heteronuclear photoassociation can proceed. In the 
language of second quantization, annihilation of an atom 
(molecule) from the i-th atomic (molecular) condensate 
is represented by the operator a i Q = a, (bo = b). The 
simplest theory has only these three levels/modes, and 
dissociation of the molecules can only take place back into 
the original atomic condensate level. To be more precise, 
molecular dissociation to noncondensate levels should be 
included, as per Fig. [TJ^b) . This situation arises because 
a condensate molecule need not dissociate back to the 
atomic condensate, but may just as well create a pair 
of atoms with equal-and-opposite momentum, since only 
total momentum is conserved. So-called rogue [HJ [21] , 
or unwanted [22], dissociation to noncondensate modes 
therefore introduces the operators a it ±^. 

The familiar mean-field equations are derived from 
the Heisenberg equation of motion for a given opera- 
tor, ihx = [x,H], which is derived from the Hamilto- 
nian (not shown), and the operators are then declared 
c-numbcrs. In a minimalist model, x represents either 
the i-th atomic amplitude a*, the molecular amplitude 
6, or the so-called anomalous density, = ai.-kC^k, 
which arises from rogue dissociation to noncondensate 
atom pairs of equal-and-opposite momentum. The sum- 
mation over k implicit to the Hamiltonian is then con- 
verted to an integral over frequency, introducing the fre- 
quency uj p — hp 2 / 3 /2m r with m r — 17111712/(1711+1112) the 
reduced mass of the atom pair. All told, the mean-field 
equations of motion are 



ia 2 



■\Slalb, 
■\Qa\b, 



(la) 
(lb) 



ib = (So-iT /2)b-ldde^f(e)A(e), (lc) 
iA(e) = eA(e)-nf(e)b, (Id) 



(a) 




la t , a 2 > 



(b) 



(c) 



6 o \ Q i <y+ 



la.a-> • 




la 1; a^> 



FIG. 1: (Color online) Few-level illustration for heteronuclear 
photoassociation. (a) Basic free-bound transition from the 
joint condensate to the electronically-excited molecular state, 
where So is the laser detuning from resonance and Q is the 
atom-molecule coupling, (b) Quasicontinuum levels account 
for photodissociation to atom pairs with equal and opposite 
momentum, where = fi/k- (c) For steady-state photodis- 
sociation, the system is effectively a two-level system with a 
lightshifted detuning S. 



where is the atom-molecule coupling, Sq is the detun- 
ing of the photoassociation laser from the |ai,a 2 ) O |e) 
transition, and the natural linewidth of the electronically 
excited state is r . A microscopic fraction of molecules 
means that Bose enhancement of spontaneous molecular 
decay [53] is negligible. Photodissociation is parameter- 
ized in terms of the kinetic energy of the atom pair he, 
the coupling 47r£ = Q/ujp^ 2 , and the continuum shape 
/(e). Whereas the standard shape for the continuum 
in quantum optics-say, autoionization theory |24)-is a 
Lorentzian, here we use a Gaussian for the sake of nu- 
merical convergence, where the width is set by the semi- 
classical size of the molecular state. 

In numerical experiments, the rate is determined from 
the time r required for the atomic probability to decay 
to ~ 1/e, 
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where Km is the rate per unit density, i.e., the rate co- 
efficient. The factor of 2 arises since we assume that 
losses are observed by monitoring one condensate or the 
other, whereas p is the total density. Note that the light- 
shift correction is obtained as follows: the detuning is set 
to an initial nonzero guess, and the time r required for 
P a ~ 1/e is determined; the detuning is updated, and 
then r is updated; repeating until r is minimized. 

We also develop a simple analytical result for the atom- 
molecule conversion rate, similar to the approach used in 
Refs. [T7] [TS] . In short, we use the adiabatic approxima- 
tion A = first to eliminate A, and again in the form 
b = to eliminate b. We thus arrive at the usual all-atom 
equations of motion with a tunable interaction 
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Which defines the (complex) resonant scattering length 
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Here the detuning and damping are now 5 — 8q + o~o 
and T = r + 7. The lightshift is er = and the 

photodissociation rate is 7 = 9(E), where 
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The rate for two-body losses from the condensate due to 
heteronuclear photoassociation is 
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Results.-Foi concreteness, we focus on photoassocia- 
tion of a joint Li-Na condensate. The atom-molecule cou- 
pling is borrowed from photoassociation of 7 Li alone [17] . 
so that fl — Q y/p/ p Q with Q = 290 x 2tt kHz and 
po = 4 x 10 12 cm 3 . Also from photoassociation of 7 Li 
alone [17], the semi-classical size of the molecular state, 
which sets the width of the Gaussian f(e), is taken to 
be L = 133cto with ciq the Bohr radius, and the natural 
linewidth of the excited state is r = 12 x 2ir MHz. 

Results are shown in Fig. [2] for p = 10 14 cm~ 3 
and 10 12 



[panel (a,c) 



cm 3 [panel (b,d)]. For lasers 



tuned to the usual lightshifted resonance [55], S = —cr , 
the numerical rate maximizes at a critical atom-molecule 
coupling f2 c (Fig. [2j dashed green lines), at a value 
that increases with decreasing density, i? max ~ 4w p for 
p = 10 14 cm" 3 and R max ~ 5w p for p = 10 12 cm~ 3 . This 
density dependence is similar to Ref. which includes 
the lightshift <tq. Correcting the lightshift according to 
a = <jq + 1.2(£!/f2 c ) 2 ro leads to saturation instead of 
maximization (Fig. [2j dot-dashed red lines) . The precip- 
itous drop in the numerical rate in panel (a) arises be- 
cause, absent the lightshift, the strongly-coupled system 
undergoes coherent oscillations between atomic conden- 
sate and photodissociated pairs, and the atomic proba- 
bility therefore takes longer to drop to P a ~ 1/e. Such 
oscillations are absent (on the timescale r) once the light- 
shift is corrected. Moreover, the corrected numerical rate 
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FIG. 2: (Color online) Many-body rate limit for heteronu- 
clear photoassociation of a joint Li-Na condensate. (a,b) The 
solid blue line is the analytical rate, and the dashed green 
(dot-dashed red) line is the numerical rate without (with) 
the unanticipated lightshift. (c,d) Photoassociation rate con- 
stant, K = R/p, where the solid blue (dot-dashed red) line is 
again the analytical (lightshifted numerical) result. The open 
red circle in panel (d) is given by K a = K m ^/(p h / p l ) 1/ ' i , 
where -Km ax is the rate constant for pi — 10 12 cm 3 and 
p h = 10 14 cm 3 . 



saturates at a value that increases slightly for increasing 
density, i? max /w p = 8.2 (8.8) for p = 10 14 ( 12 ' cm" 3 , in 
contrast to the unshifted result and Ref. [11] . The rate 
constant scales with density as K cx p -1 ^ 3 , as indicated 
by the open red circle in Fig. |2jb). Also, the analytical 
results (solid blue lines) agree best with the numerical re- 
sults for high density. Finally, R/lo p depends very weakly 
on density, so that inhomogeneous effects, which we have 
neglected, should be marginal. 

We also compare the many-body rate to the two-body 
unitary limit, which is set by the de Broglie wavelength 
of the condensate atoms according to 
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In the zero-temperature limit, the de Broglie wavelength 
is given by the mean Thomas-Fermi radius of the con- 
densate 
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where lu 3 = uj x oj v uj z defines the mean trapping frequency, 
and we assume that the frequencies can be adjusted so 
that the two condensates have identical sizes. The con- 
densate chemical potential is 
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where N is the number of atoms in the condensate, a is 
the s-wave scattering length [27] , and Lho = \jTi/{mLu) is 
the harmonic oscillator length scale. Results are shown 
in Fig. [3] for trap frequencies satisfying \J 2u> x — oj y = 
Wj/%/2 and JYj = 10 6 atoms. For loose traps, i.e., low 
density, the two-body unitary limit set by the condensate 
size is over an order of magnitude larger than the many- 
body rate limit. 

Conclusion-We have investigated the rate limit on 
photoassociation of a Bose-Einstein condensate, using an 
upgraded many-body model to obtain both numerical 
and analytical results. An unanticipated lightshift leads 
to a maximum in the numerical photoassociation loss rate 
for strong atom-molecule coupling, similar to the many- 
body model without spontaneous decay pj]. Once this 
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FIG. 3: Unitary-limited photoassociation rate constant for a 
Bose-Einstein condensate vs. trap frequency, as per the text. 
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lightshift is taken into account, the numerical rate satu- 
rates at about 9(hp 2 / 3 /2m r ), so that the many-body rate 
constant for resonant photoassociation scales with den- 
sity as K oc p -1 / 3 . This limit holds for over two decades 
in densities, whereas recent analytical results [5T] pre- 
dict this density dependence for limited densities. It also 
agrees with the rate limit obtained for combined pho- 
toassociation and Feshbach resonances [18]. Our own 
analytical approximation yields a rate constant that is 
independent of density, and agrees best with the numeri- 
cal result for dense condensates. Finally, the many-body 
rate limit has been shown explicitly to be generally more 
strict than a two-body unitary rate limit. 
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